Direct Measurement of the Photon Statistics of a Triggered Single Photon Source 
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We studied intensity fluctuations of a single photon source relying on the pulsed excitation of the 
fluorescence of a single molecule at room temperature. We directly measured the Mandel parameter 
Q(T) over 4 orders of magnitude of observation timescale T, by recording every photocount. On 
timescale of a few excitation periods, subpoissonian statistics is clearly observed and the probablility 
of two-photons events is 10 times smaller than Poissonian pulses. On longer times, blinking in the 
fluorescence, due to the molecular triplet state, produces an excess of noise. 

PACS numbers: 42.50.Dv, 03.67.Dd, 33.50.-j 



Over the past few years, there has been a growing in- 
terest for generating a regular stream of single photons 
on demand. This was mainly motivated by applications 
in the field of quantum cryptography ffl . An ideal single 
photon source (SPS) should produce light pulses contain- 
ing exactly one photon per pulse, triggered with a repe- 
tition period T rop , and delivered at the place of interest 
with 100% efficiency. For any given measurement time T, 
this source would emit exactly N = T/r rcp photons, so 
that the standard deviation AN = ^/{N 2 ) T - (A)f, = 
(( }t has to be understood as a mean value over a set 
of measurements lasting T). Such a source would then 
be virtually free of intensity fluctuations, therefore cor- 
responding to perfect intensity squeezing Q . 

A first category of SPSs already realized consists of 
sources operating at cryogenic temperature. They rely 
on optically ||, f|, ||, ^| or electrically pumped semi- 
conductor nanostructures or on the fluorescence of a two 
level system coherently prepared in its excited state [||. 
A one-atom micromaser has also been used to prepare 
arbitrary photon number states on demand ||] . However 
the collection efficiency of photons is barely higher than 
a few 10 -3 in these experiments. Due to this very strong 
attenuation, the intensity statistics are very close to a 
Poisson law at the place where the stream of photons is 
available. Another route is to realize SPSs at room tem- 
perature. In this case higher collection efficiency (around 
5%) is achieved. The existing room-temperature SPSs 
rely on the pulse saturated emission of a single 4-levels 
emitter 0. 

When the pulse duration r p is much shorter than the 
dipole radiative lifetime r ra d, such a single emitter can 
only emit one photon per pulse. This temporal control 
of the dipole excitation allows therefore to easily produce 
individual photons on demand fll2| , O ]. However, in pre- 
vious SPS realizations, little attention has been paid to 
analyse their intensity fluctuations. To address this prob- 
lem we realized a room-temperature SPS relying on the 
pulsed saturation of a single molecule embedded in a thin 
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The samples are made of cyanine dye DiICis(3) 
molecules dispersed at a concentration of about one 
molecule per 10 /im 2 into a 30 nm thick PMMA layer, 
spincoated over a microscope coverplate. The fluores- 
cence from the single molecule is excited and collected 
by the standard technique of scanning confocal optical 
microscopy Jl5fl , The molecules are non-resonantly ex- 
cited at 532 nm, with femtosecond pulses generated by 
a TkSapphire laser and frequency doubled by single pass 
propagation into a LHO3 crystal. The repetition rate, ini- 
tially at 82 MHz, is divided by a pulsepicker. The energy 
per pulse E p is adjustable by an electro-optic modulator. 
The pulse duration is about r p w 100 fs. The excitation 
light is reflected by the dichroic mirror of an inverted 
microscope, and then focused by an oil-immersion ob- 
jective (x60, NA=1.4), to form a spot of « 1/xm 2 surface 
area. The fluorescence light from the sample is collected 
by the same objective and then focused into a 30 fira di- 
ameter pinhole. After recollimation, a holographic notch 
filter removes the residual pump light. A standard Han- 
bury Brown and Twiss (HBT) setup is then used to 
split the beam and detect single photons on two iden- 
tical avalanche photodiodes. Glass filters are placed onto 
each arm to suppress parasitic crosstalk [|l6| between the 
two photodiodes. 

In order to rapidly identify single molecule emission, 
we first measure the intensity autocorrelation function 
(I(t)I(t + t)) of the fluorescence light by the standard 
Start-Stop technique with a time-to-amplitude con- 
verter ||. When a single emitter is addressed, there is 
virtually no event registered at r = 0, since a single pho- 
ton cannot be simultaneously detected on both sides of a 
beamsplitter [jlTj . The histogram shows a peak pattern at 
the pulse repetition period r rcp . As explained in Ref.[|), 
the peaks' areas allow one to infer the probabilities Ps(n) 
for the source (S), of giving n = 0,1,2 photocounts per 
excitation pulse, where 2 photons counts are due to de- 
viation from the ideal SPS. Nevertheless, this technique 
can hardly be used to extract the intensity fluctuations on 
timescale longer than a single pulse. We have therefore 



2 



chosen to record each photodetection event with a two- 
channel Time Interval Analyser computer board (Guide- 
Tech, Model GT653). Since each detection channel has 
a deadtime of 250 ns the excitation repetition rate was 
chosen to be 2 MHz. In a typical experiment, we first 
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FIG. 1: Photon counting rate R vs. energy per pulse, for a 
single molecule. The inset shows the excitation ramp E p (t), 
with _Bp lax = 13 pJ in this case. The record of the saturation 
curve was limited by the photobleaching of the dye. The 
dashed curve is a fit of the raw data according to eq.([j]) and 
the solid line is a fit after correction of triplet state excursion. 
The right scale shows the excited state population a. 

raster scan the sample at low excitation energy per pulse 
(0.5 pJ). When a single molecule is located, we apply the 
excitation energy ramp E p (t) shown on inset of Fig .[I], and 
simulteaneously record the fluorescence counts on a 1 ms 
integration time. Fig.|l| displays the fluorescence counting 
rate R vs. E p . The large intensity fluctuations are due to 
triplet state excursion of the molecule (see Fig j^) . If this 
state is not taken into account, the molecular energy lev- 
els can be modelised by a 2-level system assuming a very 
fast non-radiative relaxation between the two higher and 
the two lower energy states. The excited state population 
a at the time r p after the pulse arrival is then 

a= E v' E ^ d-e-^i^)) (1) 

(i + ^/£U)V / 

where r ra( j ~ 2.8 ns for the cyanine considered. The data 
R(Ep) are fitted by the function R = Ro x a in a two 
steps procedure. After a first fit of the raw data, all 
the points below this fit, which are attributed to triplet 
state excursion, are removed. The fit of the remaining 
set of data yields Rq = 160 x 10 3 counts/s and E sat — 
5.6 x 10~ 5 pJ. 

In order to optimize the number of emitted photons 
and avoid rapid photobleaching, we then set £™ ax to 
5.6 pJ. Such a value would correspond to a = 97%, for 
the molecule studied in Fig.|l|. During the constant max- 
imum pumping energy period of the excitation ramp, 
10 4 detection events are typically recorded before pho- 
tobleaching. Thanks to the high stability of the period 
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TABLE I: Photocount probabilities Px(n), n = 1,2, of our 
SPS (X=S), of an experimental reference source (R) and of a 
theoretical coherent source (C), the photocount statistics of 
which is affected by the detection setup, nx are the mean 
number of detected photons per pulse and V, the normalized 
variances. Negative values of V for the coherent and reference 
sources are due to a dead time effect. 

of the pulsed laser, this set of times can be synchronized 
on an excitation timebase. We then build the table of 
the number of photocounts rii = 0, 1,2 for each excita- 
tion pulse i. Photons which are delayed by more than 
10 x r ra d are considered to come from the dark counts of 
the two photodiodes, and are therefore rejected. 

The data considered hereafter corresponds to a molec- 
ular source (S) which survived during 319769 periods 
(about 160 ms) yielding 14928 recorded photons includ- 
ing 14896 single photon events, 16 two-photons events. 
We deduced P s (l) = 0.0466 and P s (2) = 5.0 x 10" 5 and 
a mean number of detected photon per pulse ns = 0.0467 
(see Table Q). The real source is considered as the su- 
perposition of an attenuated ideal SPS with an overall 
quantum efficiency 77, and a coherent source simulating 
the background, which adds a mean number of detected 
photon per pulse 7. From the measured values of Ps(l) 
and P s (2), we infer r\ « 0.0445 and 7 « 2.2 x 10~ 3 . This 
leads to a signal-to-background ratio of about 20. 

We also compared experimentally our SPS to a refer- 
ence source (R) made of attenuated pump laser pulses, 
with approximately the same mean number of detected 
photons per pulse. Quantitative tests of this reference 
source and of the detection setup are however neces- 
sary. A particular care has to be taken to the bias of 
photocount statistics due to the detection deadtimes on 
both channels. For a coherent state of light (C) con- 
taining a photon per excitation pulse, one can read- 
ily calculate the counting probability distribution and 
show that P c (0) = e- Q ,P c (l) = 2e~ a/2 (l - e~ a/2 ), 
P c (2) = (1 - e~ a/2 ) 2 and n c = 2(1 - e""/ 2 ), where 
nc is the mean number of detected photons per pulse. 
For the reference source (R), we measured jTr = 0.0462, 
Pr(1) = 0.0452, P R (2) = 50xl0~ 5 , whereas one predicts, 
for n c = n R , P c (l) = 0.0451 and P c (2) = 53 x 10~ 5 . 
The measured values are in good agreement with the pre- 
dictions, which proves that the faint TirSapphire pulses 
make a good calibration source for Poisson statistics. We 
then infer the ratio P S (2)/P R (2) = 0.10, which tells that 
the number of two photons pulses in our SPS, is 10 times 
smaller than in the reference poissonian source (R). 

In a first attempt to estimate the fluorescence intcn- 
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sity fluctuations per pulse, we considered samples of the 
data made of W excitation cycles. We introduced a nor- 
malized variance V\y defined, on the sample, by Vw = 
((An) 2 ) w /(n) w , with <(An) 2 ) w ee YZM - Ww) 2 /^, 
where n, is the number of detected photons for the pulse 
% and (n) w is the mean number of detected photon per 
pulse in the sample considered. In the very few samples 
for which (n) w = 0, Vw is not defined and is set to 1. For 
a Poisson distribution of photocounts Vw — 1, whereas 
Vw < 1 for subpoissonian distribution. 
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FIG. 2: (a) Mean number (n) w of photons detected in suc- 
cessive samples of the data of constant size W = 100 cycles 
(bin period of 50 fis). (b) Corresponding time trace of the 
normalised variance Vv, showing reduction of noise (Vw < 1) 
during the emission period of the molecule. Inset of (a) dis- 
plays the molecular level diagram. 

In order to follow the time evolution of intensity fluc- 
tuations, we then extract from the set of {n-i} all the 
successive samples of photocounts measurements of size 
W = 100, separated by a single pulse period. FigJ|(a) 
displays the mean number of detected photons per pulse 
vs. time. We clearly see random intermittency in the 
fluorescence of the molecule, due to the presence of a 
dark triplet state T\ in the molecular energy levels dia- 
gram (see inset of FigJ|(a)). At each excitation cycle the 
molecule has a small probability to jump into this non- 
fluorescent state, where it stays for a time much longer 
than the repetition period. Fig.||(b) shows, in parallel, 
the timetrace Vw(t) of the normalized variance. During 
an emission period, Vw stays below 1, and the statistics 
of the number of the detected photons per pulse is sub- 
poissonian. On the other hand, when the molecule stops 
to emit, the background light yields Vw ~ 1- If we now 
consider the whole set of data, our measurements yields 
a single value for the variance V — 0.9555. In this in- 
tensity fluctuation analysis at the level of a single pulse, 
this value of V is also directly related to the Mandel para- 
mater H Q= ({An) 2 )/(n)-l by Q = V-1 = -0.0445. 
Let us point out that due to the photodetection dead- 
time, the triggered reference source (R) also yields a sub- 
poissonian counting statistics. More precisely, for the 



coherent source (C) giving about the same mean num- 
ber nc of photons per pulse than our SPS, one predicts 
a value Q = —t%q/2 = —0.0231. This is confirmed by 
our measurements on the reference source (see Table |) . 
Nevertheless, the fluctuations of the number of detected 
photons per pulse coming out of our SPS show a clear de- 
parture from the reference coherent source. Albeit still 
limited by the quantum efficiency 77, this direct measure- 
ment of |Q| is larger by more than one order of magnitude 
compared to previous experiments |l9|, |2(|. For such a 
solid state SPS like ours, any improvement achieved in 
the light collection efficiency would therefore yield higher 
values of this subpoissonian character. We indeed ob- 
served, in preliminary experiments, an increase of the 
collection efficiency by placing the molecule at a con- 
trolled distance of a metallic mirror. 

However, the leak in the dark triplet state induces cor- 
relations between consecutive pulses. The measurement 
of the variance V of the detected photon number per 
pulse is therefore insufficient to characterize the noise 
properties of our SPS. Whereas such a characterization 
is usually infered from the record of noise power spectra, 
our photocount measurements are performed in the time 
domain. We therefore introduce, as a new variable, the 
number N(T) — J2i=i n i °f tne detected photons dur- 
ing any period of observation T — /cr rcp . The analysis 
of the fluctuations of the variable rij can be generalized 
to the variable N(T), by using the time dependent Man- 
del parameter fl) Q(T) = ((AN) 2 ) T / (N) T - 1. We 
can also define a Mandel parameter Q S (T) for the num- 
ber of photons emitted by the source in the same pe- 
riod of time T. In the case of an ideal SPS, we have 
Q = i] x Q s [fEjf . For such a source, Q s = —1, and there- 
fore Q(T) = —i], for any value of T. Fig.|| shows that 
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FIG. 3: Direct measurements of Mandel parameter Q(T) over 
short recording time T. Dashed line shows Q(T) for the equiv- 
alent coherent source (C). Inset shows Q(T) for longer time 
of observation. The solid curve is a fit to the data by a model 
accounting for intermittency in the SPS emission. 

we did observe subpoissonian intensity fluctuations on 
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timescales from T = 1 X r rop to T = 8 x r rcp , with the 
minimum value Q(r re p) = —0.0445 achieved on a single 
pulse timescale, as explained above. When we consider 
the number of detected photons on timescales larger than 
10~ 5 s, the intensity fluctuations exhibit a superpoisson- 
nian behaviour (Q(T) > 0) as shown on inset of Fig||. 
This is a direct consequence of the bunching due to the 
triplet state |23|. We developped a simple model to ac- 
count for the intermittency of the SPS emission. In this 
model, the molecule is either available for fluorescence 
and is said to be in a ON state, or it is in its triplet OFF 
state and does not fluoresce. Let us note p, the proba- 
bility per unit of time to make a ON — > OFF transition, 
and q = 1/tt the one to make the reverse OFF — > ON 
transition, where tt is the lifetime of the triplet state. 
Note that pr rep = Vise is the intersystem crossing prob- 
ability per excitation pulse. From measured values at 
the single molecule level with DiICig(3) cyanine dye 
pr rcp « 10~ 4 <C 1 and qr rcp » 2.5 x 10~ 3 <C 1. In this 
limiting case, the Mandel parameter of the source is 
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where = (p + g)r rep = Vise + t icp /t t . The Mandel 
parameter of the detected photon counts is then Q(T) = 
r\ x Q S (T). As shown on FigJ|, our data for Q(T) are well 
fitted by eq.(§) over more than 4 orders of magnitude, 
with r\ — 0.0445 (measured) and the free parameters p 
and q. The fit yields pr Tcp « 2 x 10~ 4 and tt ~ 250 /is, 
in good agreement with Ref. [Q . 

As a conclusion, the record of every photocount time 
allows one to make a direct time domain fluctuation anal- 
ysis, as presented in this Letter. This technique can 
be straightforwardly applied to other SPS, and is also 
suited to investigate photochemical propreties at the sin- 
gle molecule level. 
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